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Abstract. Observations of six luminous blue compact 
galaxies (BCGs) and two star forming companion galaxies 
were carried out with the CIGALE scanning Fabry-Perot 
interferometer attached to the ESO 3.6m telescope on La 
Silla. The observations were made in the Ha emission line 
which is prominent in BCGs. A velocity sampling of 5 
km/s and a pixel size of 0.9 arcseconds were used. In this 
paper we present the observations and the data together 
with the velocity fields and the derived rotation curves. In 
addition we provide rough estimates of the total dynami- 
cal mass and of the ionised gas mass for each galaxy. All 
galaxies display rotation, but while the companion galax- 
ies have regular velocity fields, those of the BCGs are 
complex and appear perturbed. This is the most exten- 
sive study to date of the optical velocity fields of BCGs. 
The interpretation of these results will be presented in a 
forthcoming paper (Paper II). 

Key words: galaxies: compact - galaxies: starburst - 
galaxies: kinematics and dynamics - galaxies: evolution - 
galaxies: formation - galaxies: interactions 



1. Introduction 

A Blue Compact Galaxy (BCG) is characterised by its 
blue optical colours, an Hn region like emission line spec- 
trum (therefore sometimes also referred to as an Hn 
galaxy) and compact appearance on photographic sky sur- 
vey plates. BCGs have small to intermediate sizes (as 
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measured e.g. by R25). Optical spectroscopy of BCGs in 
general reveal high star formation rates and low chemi- 
cal abundances (Searle and Sargent 1972), L equeux et al. 
1979|, Fr ench |198C| , Kunth and Sargent [L983[ Masegosa et 



al. 1994). Moreover most BCGs are rich in neutral hydro- 



gen (Thuan and Martin 1981), a requisite for the intense 
star formation generally seen. However, the gas consump- 
tion time-scale is generally much shorter than the age of 
the universe, indicating that the high star formation rate 
must be transient. Such a galaxy is commonly referred 
to as a starburst galaxy. Thus BCGs are either genuinely 
young galaxies, or old galaxies that, for some reason, have 
resume d to f orm stars at a pr odigious rate (Searle and 
Sargent |l972| , Searle et al. |1973|) . It is believed that BCGs 
undergo short (on the order of a few times 10 Myr) star- 
bursts intervened by longer (on the order of a Gyr) passive 
periods. A cyclic scenario has been proposed and may re- 
sult f rom statistical effects (Searle et al. 1973[ Gerola et al. 
1980 ). Another possibility is that supernova driven winds 
halt star formation by expelling the gas. Later, the lost gas 
might accrete back on the galaxy and create a new star- 
burst (De kel a nd Silk |1986| , Silk et al. |1987| , and Babul 



and Rees 1992| ). Other ideas incorporate galaxy interac 



tions as the triggering mechanism behind starburst activ- 
ity (cf. e.g. Lacey et al. 1993, Sanders et al. 1988). Most 



BCGs are found outside galaxy clusters and in general 
seems to be fairly isolated, although there are indications 
that Hl-companions, sometimes without obvious optical 



counterparts, may be common (Taylor 1997) 
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These different scenarios have been quite widely de- 
bated over the years, and while there is now ample evi- 
dence that most BCGs contain old stars indicating that 
the present burst is not the first one (for references see 
Paper II), there is no consensus on the process(es) that 
trigger the bursts of star formation now evident. Most ar- 
guments have been based on photometry alone. On the 
other hand the dynamics of these systems are not well ex- 
plored, still the creation of an energetic event like a sudden 
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burst of star formation is likely to have dynamical causes 
and impacts. 

To improve our understanding of the dynamics and the 
triggering mechanisms behind the starburst activity we 
have obtained Ha velocity fields, using scanning Fabry- 
Perot (FP) interferometry, of a sample of BCGs. With a 
FP it is possible to achieve a two dimensional velocity field 
with both high spatial and spectral (velocity) resolution. 
Thus we can get a much better view of the gas motions 
as compared to long slit spectroscopy. The velocity fields 
can also be used to estimate the dynamical masses of the 
galaxies, and from the Ha intensities it is possible to esti- 
mate the mass of ionised gas. Previous integral field stud- 
ies of BCGs at high spectral resolution are rare. Thuan et 
al. (1987) used Fabry-Perot interferometry to study two 
BCGs, both fainter than ours (VII Zw 403, M B = -13.5, 
and I Zw 49, M B = -17.3). While VII Zw403 showed no 
well ordered larg e scal e motion, I Zw49 to some degree did. 
Petrosian et al. ( |l997|) studied I Zw 18 (M B = -13.9) and 
found indications of solid body rotation and asymmetric 
line profiles, which they interpreted to be caused by gas 
motions in the core of the galaxy. 

In this paper we will present the FP observations of 
six luminous BCGs. These were selected to be bright in 
Ha emission. Two of the galaxies have un-catalogued but 
confirmed star forming companion galaxies and these were 
also observed. Mos t gala xies are from the sample by Berg- 
vall and Olofsson ( 1986 ) and a newer extended version of 
it. In addition, one galaxy has been taken from the cata- 
logue by Terlevich et al. ( 199l|) . In this paper (Paper I) we 
will present the observations (Sect. 2), reductions (Sect. 3) 
and the results: the derived Ha images, velocity fields and 
continuum images (Sect. 4). In Sect. 5 we describe how 
the rotation curves (RCs) were constructed and provide 
rough mass estimates based on these. In Sect. 6 we give 
comments on the velocity fields and RCs of the individual 
target galaxies. In Sect. 7 we give a short summary of the 
results presented in this paper. Throughout this paper we 
will use a Hubble constant of Hq = 75 km/s/Mpc. In pa- 
per II (Ostlin et al. 1999) we will discuss the interpretation 
of these results and their implications on the masses and 
dynamics of the galaxies and the triggering mechanism 
behind their starbursts. 



2. Observations 

Eight galaxies were selected (see Sect. 1) for observations 
on the three dark nights allocated at the ESO 3.6m tele- 
scope on La Silla, from August 30 to September 2 1995. 
Observing conditions were photometric all nights. The see- 
ing (measured at the telescope) was slightly above one arc- 
second. The exposure times ranged between 24 minutes 
(1 minute per channel) and 160 minutes (almost 7 min- 
utes per channel), plus calibration exposures. The used 
instrument was CIGALE, attached to the Cassegrain fo- 
cus of the telescope. CIGALE is basically composed of a 



focal reducer (bringing the original f/8 focal ratio of the 
Cassegrain focus to f/3), a scanning Fabry-Perot interfer- 
ometer, an interference filter (to isolate the emission line 
and suppress the sky brightness) and an IPCS detector (2- 
D photon-counting system). The basic principles of this 
instrument were described in Amram et al. (1991). The 
pixel size, projected on the sky, is 0.91" with a resulting 
field of view of ~ 4'. The IPCS, with a time resolution 
of 1/50 seconds and zero readout noise makes it possi- 
ble to scan the interferometer rapidly, avoiding problems 
with varying sky transparency, airmass and seeing dur- 
ing long exposures; and thus has several advantages over 
a CCD for this application. Narrow band interference fil- 
ters were used to isolate the region around the redshiftcd 
Ha line. Different interference filters were used for ob- 
serving galaxies with different radial velocities, see Table 
|. Only in two cases (ESO 350-IG38 and Tololo 0341-407) 
the [Nn]A6548A line might be partly transmitted. This 
line is however observed to be very weak as compared to 
Ha in the two galaxies (Bergvall & Ostlin 1999| , Terlevich 
et al. 1991) so there is no risk that this affect our results. 

In order to save observing time and to increase the 
the spectral resolution, with respect to the small velocity 
range of the observed galaxies, we decided to use a high or- 
der Fabry-Perot interferometer (hereafter simply referred 
as to FP2 for Fabry-Perot interferometer number 2, see 
Table giving at Hao a mean finesse of 10, a free spectral 
range (FSR) of 117 km/s, and a scanning step of 4.8 km/s. 
The FSR is the separation between two consecutive inter- 
ference orders and is therefore the maximum wavelength 
(velocity) range that can be observed before wavelength 
overlap occurs. The finesse is effectively the ratio between 
the FSR and the width of the Airy function. The Airy 
function (or apparatus function) is the instrumental line 
profile. Thus by scanning the FSR in a number of steps 
that equals at least twice the finesse, adequate spectral 
sampling is obtained. The spectral resolution, R, given 
in Table 1 is calculated from the width of the apparatus 
function, but for high signal to noise (S/N) the effective 
R will be higher. 

One galaxy was also observed using a lower order 
Fabry-Perot interferometer (hereafter simply referred to 
as FP1, see Table [l]), which gave oss the possibility to 
check the consistency of our results. Calibrations were ob- 
tained by scanning a narrow neon line before and after the 
observations of each galaxy. Table 1 lists some character- 
istics of the set-up used for the observations. Table 2 gives 
the target names and some parameters for the observed 
galaxies. 

3. Reductions 

Reduction of the data cubes were performed using the 
CIGALE/ADHOC software (Boulesteix |l99l| ) . The data 
reduction pro cedure has bee n extensively described in 
Amram et al. ( 1995| , 1996 , 1997 ) and references therein. The 
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Table 1. Journal of observations 
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Scanning Step 


0.36 A (16.15 km/s) 




Spectral resolution 


R > 9000 (34 km/s) 


Fabry-Perot interferometer no. 2 (FP2) 


Interference Order 


2604 at Ha (6562.78 A) 




Free Spectral Range at Hao 


117 km/s 




Finesse at Hao 


10 




Number of Scanning Steps 


24 




Scanning Step 


0.105 A (4.80 km/s) 




Spectral resolution 


R > 26000 (12 km/s) 



Table 2. List of observed galaxies. R.A. and Decl. give the right ascension (h m s) and declination (° ' ") of the targets 
for epoch 1950. my is the apparent magnitude in V. ^hoLobs is the velocity used for selecting narrow band filter for the 
observations. "Filter" is the central wavelengths of the used interference filters as measured in laboratory and corrected 
to the observing conditions: a temperature of ~ 5 °C and a focal ratio of f/8. FWHM is the full with half maximum 
of the used narrow band filter. The transmission curve for the filters have a quite squared shape, thus the full width 
at (practically) zero transmission is not much larger than the FWHM. Unless noted otherwise, all observations were 
made with FP2. 
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ESO 400-IG43B 
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Notes/references: 

1. Position and velocity reference Bergvall and Olofsson (1986). 

2. Observation made with FP2. 

3. Observation made with FP1. 

4. From an extension of the Bergvall and Olofsson sample. Position and velocity reference from NED. 

5. ESO 400-G43 and its companion ESO 400-G43B were observed simultaneously, since their separation is only 3 arcminutes. 

6. Companion toESO 400-G43, position and velocity reference Bergvall and Jorsater (1988). 

7. Position and velocity reference Terlevich et al. (1991). 

8. Quoted magnitude is in the R band. 

9. Companion to ESO 338-IG04, position and velocity reference Bergvall et al. (1998) 
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accuracy of the zero point for the wavelength calibration 
is a fraction of a channel width (< 3 km/s) over the whole 
field. OH night sky lines passing through the filter were 
subtracted by determining the emission in the field outside 



the galaxies (Laval et al. 1987). 

Since we used a high order Fabry-Perot interferometer, 
FP2, the FSR of the interferometer was of the same order 
of magnitude as the average Ha line width (~ 100 km/s) 
of the target galaxies. This means that in general we have 
no emission line free continuum. Moreover, the wings of 
the emission lines are shifted by a FSR (+1 FSR for the 
blueshifted wing and -1 FSR for the redshifted wing) and 
superimposed on the central part of the line. However this 
does not affect the velocity of the line, as long as the ve- 
locity range of the galaxy is lower than the FSR. When 
this is not the velocity jump close to one FSR could 

easily be detected and corrected for using the continuity of 
the isovelocities. The only problem is that monochromatic 
(Ha) emission could be mistaken for continuum, which in 
effect would lower the measured line intensity and over- 
estimate the continuum level. To separate the continuum 
from the body and wings of the line, we used a method 
based on fitting theoretical line profiles to the observed 
profiles. The instrumental line width (apparatus function) 
of the interferometer has a FWHM of 10±2 km/s, where 
the range reflects the (known) variation of the FWHM 
over the field of view. Thus the observed widths of the 
profiles are dominated by the intrinsic Ha line widths 
in the target galaxies. The convolution of the apparatus 
function with a theoretical Ha profile yields something 
very close to a Gaussian. The function used to fit to the 
observed profiles was a "folded" Gaussian: its wings were 
cut (at ± one half FSR from the centre of the Gaussian), 
shifted (by ±1 FSR) and added to the main body of the 
profile. The reduced datacube provide the centre of the 
central wavelength, and the sum of the continuum and 
monochromatic (Ha) emission. Thus only two parameters 
are unknown: the width and the amplitude of the Gaus- 
sian (i.e. the width and peak intensity of the Ha line). 
These two parameters are determined by a least-squares 
fit and thereby both the continuum level and the Ha in- 
tensity are determined for each pixel. To check the validity 
of this method, one galaxy (ESO 338-1G04) was observed 
with two different interferometers having different inter- 
ference order: FP1 and FP2. Using FP1, the signal mea- 
sured along the scanning sequence is separated into two 
parts: an almost constant level produced by the continuum 
light passing through the narrow band interference filter, 
and a varying part produced by the Ha line emission. 
The continuum level is the mean of all channels which do 
not contain monochromatic (i.e. Ha) signal. First of all, 
we checked that without any "data cooking", the velocity 
fields have the same shape with both FPs. The agreement 
is even better if we degrade the high resolution FP2 data 
to the spectral resolution of FP1, see also Sect. 6.3. As 
a consequence of the lines having widths comparable to 



the FSR when using FP2, the continuum level and the 
true width of the lines are not known with accuracy; and 
therefore these quantities are not presented. 

A rough flux calibration of the Ha and continuum im- 
ages of the galaxies could be made using the known instru- 
mental sensitivity. For a narrow band filter having 100% 
transmission at the observed wavelength, one photo elec- 
tron corresponds to a flux of 8.32 x 10 _17 Wm _2 s _1 , with 
an uncertainty of ±20%. Correcting for the transmission 
of the narrow band filter and the atmospheric absorption 
we could then determine the total (line plus continuum) 
flux detected. As explained above, the separation of con- 
tinuum and Ha emission is not trivial. However, in all 
cases the Ha emission dominate the detected flux. In effect 
the total Ha fluxes are known with reasonable accuracy, 
while the absolute continuum level is uncertain. Compar- 
ing with availa ble R-band photometry (e.g. Bergvall and 
Olofsson 1986 ) we conclude that the contribution from 
the continuum to the detected flux is of the order of a few 
percent, which is slightly lower than the estimate form the 
continuum fit. In Table || the derived Ha fluxes are pre- 
sented. The quoted uncertainties represent the quadratic 
sum of the intrinsic uncertainty (20 %) and the estimated 
uncertainty in determining the flux for each galaxy (10- 
30 % due to uncertainty in the continuum level and the 
exact transmission of the filter at the wavelength of the 
redshifted Ha line). The total derived Ha fluxes can be 
used to estimate A4i on , the mass of the ionised Ha emit- 
ting gas. For this we assume an electron density, n e , of 10 
cm -3 , an electron temperature of 10 000 K, a mean molec- 
ular weight, fi, of 1.23 (to include the mass of helium) and 
use the Ha effective recombination coefficient for case B 



as given by Osterbrock (1989). The assumed density is 
in the lower range of observed densities of the H II gas in 
BCGs, which usually lie in the range 10 < n c < 300 cm~ 3 



(see e.g. Bergvall 1985, and Masegosa et al. 1994) and thus 
our A^ion estimate could be regarded as an upper limit, 
perhaps one order of magnitude too high. 



4. Fabry-Perot images and velocity fields 

In Fig. 1 to 14 we present the Fabry-Perot images of 
the observed galaxies. In most cases a continuum im- 
age of each galaxy is given in the upper-left panel, and 
a monochromatic (Ha) image in the mid-left panel. The 
contour levels are usually the same in the continuum and 
monochromatic images, except that the continuum is the 
integrated intensity over the used narrow band filter while 
the monochromatic images present the total Ha flux. The 
threshold of all the images was chosen to be at 3a above 
the noise in the background. All images have linearly 
spaced contour levels. The continuum images shown have 
been smoothed with a Gaussian filter (FWHM of 3"). The 
Ha images are the full resolution maps, but in drawing the 
contours the data was slightly smoothed with a median fil- 
ter. In general, in the right panel we present the isovelocity 
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contours of the ionised gas (thick solid lines) superimposed 
on the Ha image (thin solid lines). The number of contour 
levels were here reduced to make the plots more readable. 
The heliocentric radial isovelocity lines have been drawn 
after a smoothing the original profiles with a Gaussian 
filter (FWHM 1.8" x 1.8"). The isovelocity contours are 
labelled in units of km/s. In each plot, the scale and ori- 
entation are given in the lower right corner. In general the 
relative velocities are certain to within a fraction of the 
sampling step, but there might be systematic offsets of 
a few times the FSR, since the narrow band filters used 
are four to eight FSRs broad, giving a maximum system- 
atic velocity error of a few times 100 km/s. However, for 
all galaxies there are independent velocity determinations 
available (see Table 2), which give us confidence in the 
determined systemic velocities (Table 3). 

Multicomponent decomposition 

For some of the target galaxies we see obvious double line 
profiles (or profiles with very broad and flat peaks) indi- 
cating that gas with two different velocities are present. 
When this was seen in more than a single isolated pixel 
but in a region of at least 3 by 3 pixels, and with sufficient 
S /N to rule out a noise peak, we attempted to decompose 
the velocity field into two components. This was done after 
subtracting the continuum and employing a small spectral 
smoothing with a Gaussian filter (FWHM of 15 km/s), in 
order to avoid artifacts due to noise. For this purpose a 
classical least-squares algorithm was used to fit two Gaus- 
sian profiles simultaneously to the full profile in the indi- 
vidual pixels. To reduce the number of free parameters in 
the fit, we imposed the constraint that the two Gaussians 
should have the same FWHM (as measured where they 
are well separated). The decomposition was performed in- 
dividually on single pixels, but subject to a condition of 
continuity, i.e. in determining which of the fitted Gaus- 
sian components that belong to the same dynamical large 
scale component, the velocity of either component was not 
allowed to make sudden jumps from pixel to pixel. Thus 
the velocity of one dynamical component in a certain pixel 
had to be consistent with its neighbouring pixels. 

For several galaxies, we attempted a decomposition, 
but the results are only presented when we obtained a reg- 
ular and continuous large scale secondary component. The 
failure to decompose the data cube into two components 
does not guarantee that secondary dynamical components 
are not present, only that we cannot detect them. More- 
over, we have no guarantee that the components are not 
separated by the value determined here plus the free spec- 
tral range of the interferometer. Nevertheless, the separa- 
tion between the two components are often so small than 
they should not be apparent if we had used a lower resolu- 
tion Fabry-Perot. In those cases where we have extracted 
two dynamical components, individual Ha profiles are su- 
perimposed on the velocity field of the galaxy. To make it 



more readable, one profile is usually drawn in a 2 by 2 pix- 
els square box. The abscissa of each small spectrum covers 
one FSR of 117 km/s, while the Y-amplitude is normalised 
to the peak value of the monochromatic flux in the field of 
view. In some cases where we see double or broad lines, we 
also present the line profiles of the non-decomposed veloc- 
ity field with the Y-amplitude normalised to the maximum 
intensity in the box, which makes faint parts more visible. 

For most galaxies, individual pixels or small regions 
can be found where the line profiles are double or asym- 
metric. This is an indication of small/intermediate scale 
gas motions, e.g. expanding bubbles or stellar winds. 

It cannot be excluded that several sources, overlaid on 
top of each other, have monochromatic emission within 
the narrow band filter (within a few times 100 km/s of 
the target galaxy). Then we could misinterpret a back- 
ground/foreground galaxy for a secondary dynamical com- 
ponent in the galaxy studied. This is however very unlikely 
in view of the appearance of the data. That the veloc- 
ity structure within the target galaxies (main component) 
could be a misinterpretation and rather be due to two sep- 
arate galaxies a few 100 km/s apart is extremely unlikely 
since we then would expect to see double components in 
the profiles unless there was a perfect match in velocity 
(plus a multiple of the FSR), orientation and rotational 
velocity and direction. Of course there might be more than 
two velocity components present, but we do not attempt 
a decomposition into more components, since third com- 
ponents do not display any meaningful structure over a 
scale larger than one or a few pixels. Thus third compo- 
nents (and second components in those cases where we do 
not decompose the velocity field) may reflect local small 
scale gas motions. A problem in attempting a multicom- 
ponent decomposition is that in principle any Gaussian 
can be decomposed into an infinite number of Gaussian 
components, and thus the result would be unreliable and 
model dependent. 



5. Rotation curves (RCs) 

In investigating the dynamics of the BCGs we prefer pre- 
senting a rotation curve (rather than just a position ve- 
locity cut) in order to indicate the real amplitude of the 
deprojected velocity field and to take into account data 
from a large fraction of the whole velocity field (and not 
only along a cut). However, we will see that the RCs are 
not always good tracers of the mass distribution in the 
galaxies. The RC for each galaxy has been drawn by tak- 
ing into account all velocity points within ±S degrees (in 
the sky plane) from the kinematical major axis (i.e. the 
PA). This parameter was chosen to be as big as possible, 
still allowing a regular RC. The half sector, S, used for 
each galaxy is indicated in the caption of the figure show- 
ing its RC and in Table 3. The RCs are given in Fig. 1 to 
14 at the bottom of the pages. Table 6 gives the RCs in 
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Table 3. Parameters for the rotation curves (RCs). For those galaxies where we have more than one dynamical 
component, or have drawn more than one RC, "Component" indicates which component the values refer to (see Sect. 
6) and "RC" indicates the name of the corresponding RC. The inclination (Incl) and position angle (PA) are the 
ones used in constructing (and iteratively determined from) the RCs. S is the half sector (measured from the major 
axis) of the velocity field included in calculating the RCs. Pi p h Q t and Incl p hot is the photometrically determined 
position angle and inclination respectively; the uncertainties are of the same order as for the kinematically determined 
properties. The absolute magnitudes and i?25,s are based on the quoted radial velocities in Table || and a Hubble 
parameter Hq = 75 km/s/Mpc. Whci is the determined systemic velocity in the heliocentric restframe which has a 
typical accuracy of 1 km/s. However Vhd may be systematically offset with ± one FSR (117 km/s). The i?25,s values 
are corrected for inclination (using the photometric value) and given in arcseconds and kpc. A colon after the value 
indicates that it is uncertain. 



Target Name 


Compo- 


RC v hcl 


PA 


Incl 




5' 


-P^4phot 


JncZphot 


R25,B 


M B 


note 




nent 


Fig. km/s 


O 


o 






O 


o 





"/kpc 






ESO 350-IG38 


Total 


i 




6265 


320 ± 10 


40 


± 


10 


45 


300 


35 


16.3/6.4 


-20.0 


1 


j; 


Main 


2 




6264 


320 ± 10 


40 


± 


10 


45 










2 


J3 


2:nd 


2 




6250 


140 ± 10 


40 


± 


25 


85 










3 


ESO 480-IG12 


Main 


3 




4830 


234 ± 10 


52 


± 


5 


60 


200 


64 


19.5/5.9 


-19.1 




ESO 338-IG04 


Main 


5 




2825 


60 ± 10 


55 


± 


10 


45 


55 


61 


23.5/4.2 


-18.9 


4 




Main 


I 




2816 


60 ± 10 


55 


± 


10 


55 










5 


» 


Masked 


fi 




2815 


60 ± 10 


55 


± 


10 


30 










6 




2:nd perp. 


6 


3 2824 


336 ± 5 


62 


± 


10 


20 










7 


ESO 338-IG04B 




I 




2900 


230 ± 5 


55 


± 


5 


40 


225 


58 


18.2/3.3: 


-17.5: 


8 


ESO 185-IG13 


Main 


3 




5682 


25 ± 5 


55 


± 


5 


40 


130 


42 


12.4/4.5: 


-18.7: 


8 




2:nd 


10 


5706 


225 ± 5 


50 


± 


10 


40 












ESO 400-G43 




U 


5813 


225 ± 5 


55 


± 


10 


55 


220 


40 


14.3/5.3 


-19.8 




ESO 400-G43B 






5886 


305 ± 5 


60 


± 


5 


45 


310 


73 


13.0/4.9 


-18.3 




Tololo 0341-407 


East 


Ti 


4535 


295 ± 5 


15 


± 


5 


60 






5.4/1.6: 


-16.0: 


9 




West 


il 


4560 


313 ± 5 


50 


± 


5 


50 






8.8/2.5: 


-17.0: 


9 



Notes: 

1) Non-decomposed velocity field. 

2) Primary component after decomposition of double line profiles in the centre. 

3) Secondary component after decomposition of double line profiles in the centre. 

4) Total velocity field from FP2 data. 

5) Total velocity field from FP1 data. 

6) FP2 data, parts of the velocity field have been masked, cf. Sect. 6.3. 

7) Secondary perpendicular component, which has no obvious counterpart in broad band images, cf. Sect. 6.3. FP2 data. 

8) R,2b,b and Mb estimated from V-band data assuming B — V — 0.5; this value is uncertain. 

9) Tololo 0341-407 has been decomposed into two components. Total absolute R-band magnitude: Ma = —18.2. The 
kinematical inclination and position angle have been used also for calculating -fibs.s- In addition B,25,b has been estimated 
from an R-band luminosity profile assuming B — R = 0.8. Thus the values of i?25,s for this object are uncertain. 



tabular form (this Table is only published electronically )[]. 
The "cloud" of small points seen in each RC, are all the 
velocity points within ±30° of the major axis in the sky 
plane. Although these points represent only a portion of 
the data (usually more data points are included in con- 
structing the RC since normally S > 30°), they give a 
fair idea of the amount, dispersion and quality of the data 
and the difference between the velocity on the major axis 
and the mean velocity. If the discrepancy is large, it means 
that circular motion is not likely. The error bars represent 

1 Table 6 is only available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 



http: / / cdsweb.u-strasbg.fr/Abstract.html 


, it is a 


at 


ftp : // ft p. iap . fr /pub/ from_uscrs /ostlin / 


Articles 



the ±lcr dispersion in the velocity at each radius in the 
RC, and is a combination of the intrinsic dispersion and 
observational and reduction errors. In general the intrinsic 
dispersion is larger than the observational errors. Thus the 
dispersion in the RC is in general caused by real irregular- 
ities in the velocity fields. The rotation velocity scale has 
not been adjusted by the cosmological correction (1 + z). 

The RCs assume axisymmetric objects with circular 
rotation and are sensitive to the choice of inclination 
(Incl), position angle (PA) and dynamical centre. The 
PA was determined from the orientation of the velocity 
field with a typical accuracy of 5 to 10 degrees. The incli- 
nation was determined as to minimise the residuals and 
the dispersion in the RC. The centre coordinates and sys- 
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temic velocity were chosen to give a RC with good agree- 
ment between the receding and approaching sides and to 
minimise the dispersion. The centre coordinates could be 
determined with a typical accuracy of half a pixel and in 
general the displacement between the dynamical centre 
and the broad band photometric centre (determined from 
the 1st statistical moment) is within one pixel. The incli- 
nation is in general the most uncertain parameter. Since 
many velocity fields look perturbed, a wide range in Incl 
may give comparably good fits overall. A typical accuracy 
for Incl is 10°. Secondary components have in general less 
well determined Incl. The uncertainties in Incl will enter 
in all mass estimates based on the RCs. 

The RC gives V(R), the rotational velocity as a func- 
tion of radius, and in essence the rotational velocity is the 
velocity in the sky plane divided by sin( Incl). Thus the 
lower Incl, the more sensitive will the derived velocities 
will be to uncertainties in Incl. The possible occurrence 
of warps or oval distortions and irregularities limit the 
validity of the RC, since its derivation is based on the 
assumption of a circular motions. 

A rough estimate of the dynamical mass can be made 
using the simple model by Lequeux (1983) in which the 
mass within a radius R is: 



M(R) = f x Rx V 2 (R) x G~ 



(1) 



where V(R) is the rotational velocity at R, G the 
gravitational constant and / is a constant which has a 
value between 0.5 and 1.0. This formula is valid for any 
galaxy (in equilibrium) supported by rotation. For a disc 
with flat RC / = 0.6 while for a spherical distribution 
(e.g. a galaxy dominated by a dark halo) / = 1.0. For a 
disc with Keplerian decreasing RC outside R, f = 0.5 . 
Thus, according to this model, for any rotating galaxy / 
should lie in the range: / = 0.5 to 1.0, independent of the 
presence of a massive halo. Mass estimates based on this 
equation are given in Table 4, where we assumed / = 0.8. 
The assumption / = 0.8 is not based on any physical rea- 
son, but was chosen simply to lie in the middle of the 
allowed interval. In addition to the intrinsic uncertainty 
in this model, all the uncertainties above affects the accu- 
rateness of this estimate. In Paper II we will discuss these, 
and more refined, mass estimates of the observed galax- 
ies. Of course, the mass estimate provided by Eq. (1), will 
only be a good approximation of the true dynamical mass 
if the galaxy is supported by rotation. If on the other hand 
the galaxy is mainly supported by random motions, the 
presented mass will be a severe underestimate. 

The RC is based on the assumption of circular rota- 
tion. If this is not true, the RC will not give a good de- 
scription of the dynamics and mass of a galaxy. Never- 
theless, even when the velocity field is perturbed and the 
derived RC looks weird (as for ESO 338-IG04, Fig. |), the 
RC gives some insights to the dynamics of a galaxy. Even 
the failure of constructing a symmetric and tight RC is 



interesting, since it indicates that the system is complex 
or perturbed. In Table 3 we give some parameters for the 
RCs. In some cases where we have extracted more than 
one component, we provide information for both. We also 
provide some photometric information, like Inclphot and 
PAphot, the photometric inclination and position angle, 
respectively; and i?25,Sj the radius at which the B surface 
brightness drops to 25 magnitudes per square arcsecond, 



corrected for inclination (cf. e.g. Bergvall et al. 1999) and 



Galactic reddening. This information is not used in itself 
in the present investigation but provide complementary in- 
formation. The i?25,s is given in the RCs to give an idea 
of the size of the galaxy and the extent of the RC. When 
we did not have B-band data, we scaled V- or R-band data 
assuming crudely B — V = 0.5 and V — R = 0.35. The po- 
sition angle and inclination derived from kinematics and 
photometry do not always agree which could be due to 
dynamical disturbances and instabilities in the systems 
or internal absorption. Anyway, the only thing we used 
PAphot and Inclphot for in this investigation was to calcu- 
late i?25,B- Detailed surface photometry of some galaxies 
in this sample will be presented separately (Bergvall and 
Qstlin |1999|). 



6. Notes on the individual targets 

In this section we describe the images, velocity fields and 
rotation curves of the individual targets. In paper II, the 
characteristics of the observed velocity fields and RCs will 
be analysed and interpreted. In Table 3 we list the inclina- 
tions and position angles, derived from both photometry 
and kinematics. The kinematical centres have been com- 
pared with those from photometry, based on a first mo- 
ment fit and the centre of outer isophotes. Unless noted 
otherwise, the agreement between the kinematic and pho- 
tometric centres is within one pixel. 

6.1. ESO 350-IG38 

Our data (Fig. 0, ||) and broad band images show that the 
morphology of this object is complex, apparently involv- 
ing three nuclei/hot-spots. The outer Ha and continuum 
isophotes appear fairly round. However, deep broad band 
images however reveal an asymmetric morphology at all 
radii. The kinematical centre of this galaxy lies between 
the north-west and south nuclei, approximately in the cen- 
tre of the outer Ha and continuum isophotes. The isove- 
locity contours are very squeezed in the centre indicating 
a rapid increase in the rotational velocity. At larger radii 
(1 to 2 arcseconds) a plateau is observed on each side fol- 
lowed by a velocity drop after which the RC stays flat at a 
rather low rotational velocity. Note that the steep central 
velocity gradient is not necessarily continuous, but could 
also be explained by Ha emitting blobs with different ve- 
locity tumbling around each other. The deformations of 
the isovelocity contours (strongest in the north-east re- 
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Table 4. Order of magnitude mass estimates. The first 
column gives the target name. In those cases where the 
target seems to be composed of more than one dynamical 
component, column 2 indicates which component is con- 
sidered: "Total" means that it is the total non-decomposed 
velocity field, "Main" means that it is the main component 
of the galaxy, "2:nd" that it is the secondary (normally 
weaker), and " Comp" means that it is a companion galaxy. 
Column 3 gives the rotational velocity at the radius R, at 
which the mass is calculated. Column 5 contains the mass 
estimate, M. , where we have assumed / = 0.8 in Eq. (1). 
See Sect. 5 for a discussion on the accuracy of these mass 
estimates. Sometimes there are more than one mass esti- 
mate for each galaxy (component); this is when the mass 
have been calculated at several radii, e.g. at the radius of 
maximum rotational velocity (indicated by V max ) and at 
the last point in the RC (indicated by R m ax), if these do 
not coincide. In the other cases, the maximum rotational 
velocity (V max ) occurs at the last point (R m ax) in the RC 
and these are also the values of V(R) and R given in the 
Table below. The note "Both" indicates that the mass has 
been evaluated at the last point in the RC where there is 
data both from the receding and approaching sides. 



1993). Thus the full line pro- 



files are much broader than our FSR, and the problems 
discussed in Sect. 3 are more serious for this galaxy. 

The derived RC shows a rapid rise followed by a Keple- 
rian like decline shortly after, then levelling out at a con- 
stant velocity of 30 km/s for radii larger than five arc- 



Table 5. Estimated Ha fluxes (/Ha), luminosities (Ln a ) 
and ionised gas masses(A^ion)- The Ha luminosity was 
calculated using the quoted radial velocities in Table |[ 
The mass was calculated assuming an electron density of 
n c = 10 cm -3 , an electron temperature T c = 10 000K, 
and a mean molecular weight of fj, = 1.23. The Ha re- 
combination coefficient was taken from Osterbrock, Case 



B (1989) 



Target 


/Ho 

10" 17 W/m 2 


Lhc 

10 33 W 


Mion 

1O 6 X 


ESO 350-IG38 


350 ± 180 


280 


900 


ESO 480-IG12 


130 ± 50 


62 


200 


ESO 338-IG04 


290 ± 60 


49 


160 


ESO 338-IG04B 


24 ± 5 


4 


14 


ESO 185-G13 


120 ± 35 


81 


260 


ESO 400-G43 


240 ± 70 


170 


550 


ESO 400-G43B 


28 ± 6 


20 


65 


Tololo 0341-407 E 


45 ± 13 


19 


61 


Tololo 0341-407 W 


33 ± 9 


14 


45 



Target 


Compo- 


V(R) 


R 


M 


note 




nent 


km/s 


kpc 


10 6 X Q 




ESO 350-IG38 


Total 


95 


0.2 


340 


Vmax 


JJ 


Total 


30 


5.4 


900 


Rmax 


ESO 350-IG38 


Main 


41 


4.7 


1470 




ESO 350-IG38 


2:nd 


85 


1.0 


1340 




ESO 480-IG12 


Main 


120 


2.8 


7480 


Both 


jj 


Main 


140 


4.5 


16370 


Vmax 


jj 


Main 


124 


5.8 


16550 


Rmax 


ESO 338-IG04 


Masked 


49 


2.0 


890 


Vmax 


ESO 338-IG04 


Masked 


35 


2.5 


570 


Rmax 


ESO 338-IG04 


2:nd 


19 


0.8 


54 


Vmax 


JJ 




12 


1.9 


51 


Rmax 


ESO 338-IG04B 


Comp. 


91 


3.5 


5380 




ESO 185-G13 


Main 


51 


3.8 


1830 




ESO 185-G13 


2:nd 


26 


2.7 


340 




ESO 400-G43 




54 


1.1 


610 


Vmax 


JJ 




25 


5.6 


650 


Both 


JJ 




12 


10.2 


270 


Rmax 


ESO 400-G43B 


Comp. 


62 


4.3 


3070 




Tololo 0341-407 


East 


24 


1.6 


170 




Tololo 0341-407 


West 


27 


1.9 


260 





gion) could indicate the presence of a spiral arm or a warp. 
The emission lines of this galaxy are broad, on the order 
of 200 km/s (Vader et al. ~ 



seconds. Between approximately one and two kpc the RC 
falls slightly faster than the Keplerian prediction, but this 
is within the uncertainties. The receding and approaching 
sides are in reasonable agreement. 

A close look a the line profiles in the centre revealed 
double lines which made us attempt a decomposition of 
the velocity field (see Fig |2|). This yielded a regular sec- 
ondary component in the central part only. After decom- 
posing the velocity field the RC looks very different (see 
Fig. ||). The central high velocity part in Figgis replaced 
by a secondary counter rotating component spinning at 
high velocity. To get good agreement between the reced- 
ing and approaching sides of the secondary component, 
its systemic velocity was adjusted. Note that the RC of 
the secondary component is based on a total of nine pix- 
els only. The centre of the primary component nearly co- 
incides with that of the total non-decomposed velocity 
field. The primary component has a nearly constant RC 
with a velocity of w 40 km/s. This differs from the non- 
decomposed value of 30 km/s because it is based on fitted, 
i.e smoothed, data, and slightly different (1 pixel) centre 
coordinates. Moreover, since fitting Gaussian profiles to 
the data requires higher S/N, the radial extent of this RC 
is somewhat smaller than for the non-decomposed velocity 
field. A word of caution is necessary here: the velocity dif- 
ference between the primary and secondary components is 
close to the FSR of the used FP2, meaning that the lines 
almost overlap. In effect the relative velocity of the sec- 
ondary component with respect to the primary, is uncer- 
tain. Therefore we cannot exclude that the counter rotat- 
ing component is an artifact. However there is a clear sig- 
nature of asymmetric line profiles in the centre and some 
sort of multicomponent gas is needed. Observations with a 
FP with higher FSR (e.g. FP1) should reveal if the second 
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component represents a counter rotating disc or not. At 
this stage we urge the reader to view the RC of the sec- 
ond component in Fig. || as one possible interpretation. 
If our interpretation is correct, the mass of the secondary 
component is on the order of 1O 9 .M0. 

This galaxy has the greatest Ha luminosity in the sam- 
ple and the ionised gas mass may be as great as 10 9 .M©, 
comparable to the dynamical mass estimate in Table 4. 

6.2. ESO 480-IG12 

This galaxy also has an overall irregular morphology. It is 
located in front of what appears to be a background cluster 
of galaxies. Three of the galaxies closest to ESO 480TG12 
have measured redshifts that are much higher than that of 
ESO 480-IG12. Furthermore, none of these objects show 
any Ha emission at the rcdshift of ESO 480-IG12. Surpris- 
ingly though ESO 480-IG12 has faint warp- like extensions 
apparently aligned with a chain of background galaxies. 

The Ha emission is concentrated to two, or perhaps 
three, central Hn regions. In between the Hn regions the 
isovelocity contours are squeezed. Double and very broad 
components are observed in the emission line profiles, es- 
pecially west of the centre (see Fig. || and ||). Two com- 
ponents have been extracted with quite different inten- 
sity levels. If we do not decompose the profiles and sim- 
ply compute the velocity from the total Ha line, the first 
component clearly dominates in the brightest part, but in 
fainter regions west of the centre, unreasonable features 
appear like the confusing situation with three different 
symmetry axes. Anyway, this would not affect the RC of 
the primary component since the regions with double lines 
are close to the minor axis. For the primary component, 
the velocity field shows a strong gradient along the major 
axis; but with a plateau in the north-west region. This 
coincides with the western Hn region and roughly with 
the region where the double features in the lines are most 
pronounced, and thus close to where the second compo- 
nent has its maximum intensity. The kinematical centre 
is well defined along the major axis and coincides roughly 
with the continuum peak intensity. Due to the plateau, 
the kinematical centre is less well determined along the 
minor axis, however different choices of centre along this 
axis gives consistent RCs. In essence the RC does not sen- 
sitively depend on the centre coordinates or the decompo- 
sition model. 

The RC derived for the main component shows solid 
body rotation out to a radius of « 3 kpc, after which 
it levels out, although it is here based on the receding 
side only. The RC has a rather large dispersion indicating 
that the physical situation might be more complex than 
the assumed disk anatomy. For the second component, we 
could not obtain any sensible RC. However, the duality 
of the velocity peaks in the western part of the galaxy 
are very significant, and some sort of multicomponent gas 
is needed to explain the data. The ionised gas mass is 



of the order lO 8 ^©, thus significantly smaller than the 
dynamical mass estimate which is of the order lO lo yVf . 



6.3. ESO 338-IG04 

This galaxy is also well known as Tololo 1924-416 or SCHG 
1924-416. Its photometric properties have been quite ex- 



tensively discussed by Bergvall (1985) and Ostlin et al. 
(1998). This galaxy has been observed with two different 
Fabry-Perot interferometers (FP1 and FP2, see Table 1) 
and thus we had the opportunity to check the consistency 
of the analysis (cf. Fig. || and |?|). First, we successfully 
checked that the flux in the field star superimposed on 
the western part of the galaxy and the flux in the Ha 
emission line regions were consistent between both ob- 
servations. Secondly, we checked that the superimposed 
bright field star does not seriously affect the velocity field. 
Thirdly, we confirm that we find the same shape for both 
velocity fields. 

The Ha emission is concentrated to an extended bright 
central starburst region. In addition there is diffuse Ha 
emission extending in a tail towards the east. There are 
also suggestions of a small Ha arm emanating towards the 
south from the western side of the starburst. The velocity 
field is very irregular and does not contain a single axis of 
symmetry. Moreover the velocity gradient is steep in the 
western parts and roughly east-west orientated; while in 
the eastern regions the gradient is much lower, not always 
positive and lacking a well defined position angle. The 
eastern extended tail has almost no velocity gradient. 

East of the centre, just at the border of the starburst 
region, we observe what appears to be the superposition 
of two patterns with different orientation (see Fig. ||). We 
will refer to the hypothetical component east of the centre 
as the perpendicular component, since its PA is roughly 
perpendicular to the major axis of the galaxy. However, 
we do not observe any double component in the profiles 
(in either FP2 or FP1 data). Anyway, we tried to decom- 
pose the velocity field without any conclusive results. This 
could mean that if there are two components present, the 
lincwidth is too large with respect to their velocity sep- 
aration or/and where the components overlap, one only 
sees the one with the strongest Ha emission. 

In Fig. |7] we show a map of the velocity dispersion as 
derived from the FWHM linewidth of the FP1 data. The 
velocity dispersion has a fairly constant level of « 100 
km/s. Where the major axis of the perpendicular compo- 
nent cross the major axis for the whole galaxy the veloc- 
ity dispersion is higher and peaks at 160 km/s. Still, the 
shape of the Ha line is consistent with one single broad 
component. 

It is not unambiguous how to derive a RC for this 
galaxy. The different RCs are however consistent in the 
way that they all are very irregular, signifying a non equi- 
librium system. The general feature is an approaching side 
with continuous steeply rising velocity, and a receding side 
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with very small velocity gradient. This means that the as- 
sumption of a regularly rotating disc in equilibrium must 
be far from reality in this galaxy. In effect the kinematical 
centre is not well defined. 

As a cure we tried to mask away certain points of the 
velocity field (those in the eastern arm and the approach- 
ing side of the perpendicular component) to check if we 
could obtain a more regular RC. The result is shown in 
Fig. ||a, but the RC is still far from regular. Moreover, the 
decline of the average RC outside 2 kpc is faster than the 
Keplerian case, which is not necessary significant since the 
velocity field is so perturbed. 

For the hypothetical perpendicular component we 
could extract a regular RC (Fig. gb) only when using a 
narrow sector (S < 30°). The implied rotational velocity 
is very small, but the regularity and the agreement be- 
tween the approaching and receding sides still makes this 
component well defined. 

The RC obtained from the FP1 data (Fig. 0) largely 
agrees with that from FP2 data (Fig. ||), but the rota- 
tional velocity of the receding side is » 10 km/s larger for 
the FP1 data. The differences can be attributed to three 
effects: Firstly, the different spectral resolution of FP1 and 
FP2. Secondly, the slightly different choice of centre coor- 
dinates (an independent determination was made for the 
FP1 data). Thirdly, the FP1 data is somewhat deeper and 
reach larger radii, and since the velocity field is depro- 
jected these outer points enter also at smaller radii. In 
view of this, the agreement is very good. 

A last note on this galaxy is that here it is question- 
able that the expression " rotation curve" is adequate since 
there appears to be no regular rotation present. Still, the 
presented RCs provide important information on the com- 
plexity of the system. However, the dynamical mass esti- 
mate in Table 4 should be viewed sceptically. The ionised 
gas mass is of the order 1O 8 A40. 

6.1 ESO 338-IG04 B 

This galaxy is a physical companion of ESO 338-IG04 and 
is located approximately six arcminutes south-west of it, 
corresponding to a projected distance of w 70 kpc. The 
Ha image contains at least four bright Hn regions (Fig. 
|8|). The galaxy has an irregular, somewhat "bent" shape 
in Ha but the velocity field is still quite regular. Along 
the north-west side there are twists in the isovelocity con- 
tours which may indicate the presence of a spiral arm. 
Just north-east of the brightest H n region there is a steep 
velocity gradient which we interpret as the kinematical 
centre, and which coincides with the centre of broad band 
images. The RC shows good agreement between the two 
sides and a gradual flattening with increasing radius. The 
estimated dynamical mass is 5 • 10 9 A^© and the ionised 
gas mass is of the order of 1Q 7 Mq. 

On broad band CCD images, about 40" north-west of 
the target galaxy, there is a small low surface brightness 



galaxy, which is not detected in the monochromatic Ha 
images. 

6.5. ESO 185-IG13 

This galaxy has a nearly round shape in the continuum 
but a more complex Ha morphology and velocity field 
(Fig. ||) . Broad band images reveal arms and a > 30" long 
tail extending to the north-east. The Ha emission is strong 
over a large area with two peaks in the surface brightness. 
The velocity field appears squeezed in the north-west re- 
gion, between the two peaks in the Ha surface brightness. 
About 1.5 arcminutes south-west of the galaxy there are 
two galaxies present, none of which is detected in Ha at 
the redshift of ESO 185-IG13. 

Double components with rather different intensity 
have been extracted from the original lines. Both com- 
ponents have similar position angles but the gradients are 
reversed, i.e. the secondary component is counter-rotating 
with respect to the primary (Fig. gand 10). The brightest 
component is very similar to the total non-decomposed 
velocity field (compare the upper left and right panels in 
Fig. g). The second component is more diffuse and regu- 
lar than the first one. While the centre of the secondary 
component coincides with the continuum peak, the kine- 
matical centre of the primary component is slightly offset 
from it (« 1" to the north). 

The RC of the primary component shows a rapid in- 
crease followed by an almost flat slowly rising part. The 
agreement between the approaching and receding sides is 
good. The secondary component also yields a regular RC, 
with a low rotational velocity though. The regularity of the 
RC of the secondary component suggests that it is caused 
by a dynamically well defined object, e.g. a counter rotat- 
ing disc. The estimated dynamical mass of ESO 185-IG13 
is ~ 2 • 10 9 Mq and the ionised gas mass is of the order of 
10 8 A4©. The dynamical mass of the secondary component 
is - 3 • 1O 8 7W 

6.6. ESO 400G-43 

This galaxy presents several bright H11 regions. Note that 
the northernmost region is not a field star. The south- 
ern H 11 region complex has a ring-like morphology. In the 
southern H 11 region complex we observe a symmetric and 
regular velocity field with velocity plateau's followed by 
decreasing velocity on both sides, characteristic of a disc 
with circular differential rotation. At larger radius (> 5") 
the north-east and south-west regions start to behave dif- 
ferently. In the south-west we observe the continuity of 
the central velocity field, although the isovelocity contours 
are somewhat boxy. The region north-east of the centre is 
more intriguing because we observe a shift of the major 
axis position angle, rotating from approximately 45° to- 
wards north, which might indicate the presence of a warp. 
An alternative interpretation is that the twisted isoveloc- 
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ity contours north of the north-east velocity plateau is due 
to the presence of a spiral arm. The isovelocity contours 
in the northern part of the velocity field seems to indicate 
the presence of local motions, and asymmetric line profiles 
are present. The kinematical centre is well defined but is 
slightly offset (« 1" towards south-west) from the peak 
continuum emission. Due to the irregular velocity field in 
the north east, the inclination is quite uncertain for this 
galaxy. 

The RC has a strange behaviour. Following the initial 
solid body rise to the plateau's, there is a rapid decline 
after which the RC levels out and stays flat out to a ra- 
dius of 15 arcseconds, after which the approaching side 
(the only one that has Ha signal) declines further. The 
remarkable thing is that after the maximum rotational 
velocity the decline is faster than for the Keplerian case. 
This is unphysical for an equilibrium disc, thus indicating 
that the velocity field is distorted. The "super-Keplerian" 
decline can be eased somewhat, but not completely, by 
allowing the position angle to be variable. In the south 
west regions at a radius greater than 5" double features 
in the Ha lines appear at low S/N. However, this feature 
is present in many consecutive pixels and is in total sig- 
nificant, which means that a secondary component might 
be present here. In the centre there are no signs of dou- 
ble lines and our attempts to decompose the velocity field 
failed. Due to the super-Keplerian behaviour of the RC 
the estimated mass apparently decreases with increasing 
radius (see Table 4). Thus the mass stated in Table 4 is 
very uncertain. The ionised gas mass is of the order of a 
few times lO 8 ^©, comparable to the dynamic estimate. 



6.7. ESO 400G-43 B 

This quite low surface brightness galaxy (Fig. |l^) is a 
physical companion of ESO 400G-43 and is located ap- 
proximately three arcminutes east /south-east of it (Berg- 
vall and Jorsater 198? ), corresponding to a projected dis- 
tance of w 70kpc. ESO 400G-43B has one extended central 
emission line region, with peak intensity in the north-west. 
The isovelocity contours shows axisymmetric twists south- 
east and (to a lesser extent) north-east of the major axis, 
perhaps indicating the presence of spiral arms. The shape 
of the RC may be affected by dust extinction, due to the 
high inclination. However, the best fitting kinematical in- 
clination is significantly lower than the photometric incli- 
nation derived from ellipse-fitting of the outer isophotes, 
the reason for which is not well understood. In the pre- 
sented RC we have applied a small correction for internal 
extinction in the central parts, see Fig. 12. The estimated 



dynamical mass is ~ 3 • 10 9 AiQ and the ionised gas mass 
is of the order of a few times 10 7 Mq. 



6.8. Tololo 0341-4 07 

Even if this object has a short observing time (see Table 
2) it was possible to extract a high S/N velocity field. Our 
data indicate that this object is in fact composed of two 
well separated galaxies: the eastern one being more lumi- 
nous in Ha than the western one (Fig. [l3| and [lj]). The ve- 
locity field of the eastern component displays a rather con- 
stant velocity indicating a close to face on orientation. The 
western component is more inclined and shows more obvi- 
ous rotation. Asymmetries in the line profiles suggested a 
decomposition of the velocity field into two components, 
one which is dominated by the eastern galaxy and the 
other by the western one 

(Fig. H and 0). If the 

connec- 
tion between the first component in the eastern and west- 
ern galaxy is physical is not certain; and the same is true 
for the secondary components. Here, the notation with 
a primary and secondary component may be somewhat 
confusing, since the secondary component of the western 
galaxy is actually the most luminous. Hereafter, when we 
mention the "first" component of the eastern or western 
galaxy we refer to the strongest (most luminous) compo- 
nent, and vice versa for the secondary components. The 
secondary component of the eastern galaxy has a constant 
velocity, which suggests that this component belongs to 
the same galaxy, although the velocity is more similar 
the western galaxy. The second component of the western 
galaxy has a major axis position angle roughly perpendic- 
ular to the one of the first component, which suggest that 
it is dynamically distinct from it. The primary compo- 
nent of the eastern galaxy has a solid body RC with quite 
high dispersion. The primary component of the western 
galaxy gives a well determined RC, perhaps levelling out 
at large radii. The kinematical centres roughly coincide 
with the peaks in the continuum emission, although this 
is hard to determine for the eastern galaxy. The secondary 
components do not produce regular RCs. The eastern and 
western galaxies both have estimated dynamical masses 
around ~ 2 • lO 8 ^©, and the ionised gas masses are of 
the order of a few times 10 7 A4q. 

7. Summary 

We have presented velocity fields for a sample of six rela- 
tively luminous blue compact galaxies, and two compan- 
ion galaxies. In general the velocity fields, in spite of high 
S/N, appear irregular and distorted, except for the two 
companion galaxies included in the sample. We have also 
shown the derived rotation curves (RCs) for the simplest 
choice of input parameters. That some RCs appear strange 
or non-uniform indicates that the simple assumption of a 
regular warp-free disc with circular rotation is not valid in 
general. Still the information in the RCs are valuable to 
get a feeling of the dynamics of these systems and to get a 
rough estimate of the mass contained within the observed 
extent of the galaxies. The estimated masses range from a 
few times 10 8 to ~ 10 10 M.q. Secondary components have 
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smaller estimated masses (when a RC could be derived). 
These mass estimates are approximate and based on the 
assumption that the galaxies are supported by rotation 
alone. We will show in Paper II that some of the galaxies 
cannot be primarily rotationally supported. 

Using the integrated Ha fluxes, the masses of ionised 
gas were estimated to lie in the range 10 7 to lO 9 ^©. This 
amounts to between less than 1% and more than 50% 
of the estimated dynamical masses. The two companion 
galaxies had the lowest fraction of ionised gas mass to 
dynamical mass, and the galaxies with the most peculiar 
RCs (ESO 350-IG38, ESO 338-IG04 and ESO 400-G43) 
apparently have the largest fraction of their mass in the 
form of ionised gas. 

This paper presents the most extensive study of the 
central dynamics and optical velocity fields of BCGs as 
yet. In addition, all the observations have been done with 
the same instrument and the sample selection, reductions 
and analysis are homogeneous. In a forthcoming paper 
(Paper II) we will discuss the interpretation of the veloc- 
ity fields and rotation curves in terms of more detailed 
mass models. In addition we will discuss what the results 
presented here tell us about the triggering mechanism for 
the starburst activity present in these galaxies. 
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Fig. 1. ESO 350-IG38. Upper left: Continuum isointcnsity image. Mid-left: Monochromatic (Ha) isointensity map. 
The (relative) step between two consecutive levels is the same in both images. Upper right: Isovelocity contours of 
the ionised gas (thick solid lines) superimposed on the Ha image (thin solid lines). The spatial scale in each image 
is shown in the lower right corner. North is up, east is to the left. Bottom: Rotation curve (RC) derived from the 
velocity field and based on S — 45°, incl = 40° and PA = 320°. The receding side is indicated by (J) and the 
approaching side by ©. The weighted average of the two sides is indicated by • with the lengths of the errorbars 
equal to the Ha dispersion of the data points within the sector. Velocity points which are based on only one or two 
data Doints arc enclosed bv brackets (). The curve (solid line) is a smoothed SDlincfit to the average RC. The cloud 
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Fig. 2. ESO 350-IG38 continued. Upper left: Line profiles in the centre of the non-decomposed velocity field plus 
isovelocity contours of the total, non-decomposed, velocity field. Each box corresponds to one pixel. The Y-axis scale 
is normalised to the peak intensity in each pixel. Mid-left: The line profiles in the centre of the two fitted components 
after decomposition of the velocity field. Each box corresponds to one pixel and the intensity scale is normalised to 
the pixel with the brightest Ha flux. Upper right: Isovelocity contours of the first component of the ionised gas 
superimposed on the Ha emission from the first component. Bottom: Rotation curve (RC) based on the decomposed 
velocity field. The usual symbols corresponds to the first component, and are based on S = 45°, incl = 40° and 
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Fig. 3. ESO 480-IG12: Upper left: Continuum isointensity image. Mid-left: First monochromatic (Ha) component 
plotted with the same step as in the continuum image. The isovelocity contours of the first component are overlaid (thick 
solid lines). Upper right: Isovelocity contours of the first component of the ionised gas (thick lines) superimposed 
on the two Gaussian components of the two velocity fields (thin lines) plotted within squares of 2 pixels. The spatial 
scale in each image is shown in the lower right corner. North is up, east is to the left. Bottom: Rotation curve (RC) 
based on S — 60°, incl = 52° and PA = 234°. For further explanation of what is shown in the RC see the caption 
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Fig. 4. ESO 480-IG12 continued. Upper left: Monochromatic (Ha) isointensity image of the total (non-decomposed) 
velocity field. Bottom left: Second monochromatic component plotted with the same step as in the continuum of the 
first component images (Fig. [|). Right: Isovelocity contours of the second component of the ionised gas (thick lines) 
superimposed on the non-decomposed line profiles (thin lines) plotted within squares of 2 x 2 pixels. The intensity 
scale (Y-axis) in each box with line profiles is normalised to the peak intensity in that box. The spatial scale in each 
image is shown in the lower right corner. North is up, east is to the left. 
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Fig. 5. ESO 338-IG04: Upper left: Continuum isointensity image. Note the bright field star superimposed on the 
western part of the galaxy. Mid-left: Monochromatic (Ha) isointensity map. The step between two consecutive levels 
is the same for the continuum and monochromatic images. Upper right: Isovelocity contours of the ionised gas (thick 
solid lines) superimposed on the monochromatic image (thin solid lines). The step between two consecutive contours 
is 10 times larger than the step used in the mid-left image. The spatial scale in each image is shown in the lower right 
corner. North is up, east is to the left. Bottom: Rotation curve (RC) based on S = 45°, incl = 55° and PA = 60°. 



18 



Ostlin et al. 







- o 

I 00 



o 

CO 







i 9 

co CM 



-+->_ 
O 

_o 

^ o 

o 
c 
o 



o 



o 







Radius (kpc) 



2 



1 1 — 

ESO 338-IG04 

Both sides 
q Approaching side 
Receding side 



CP 








~T 1 

o o o 
o 



5 10 
Radius (arcsec) 

Radius (kpc) 
1 



ES0338-IG04 




Receding side 
q Approaching side 
m Both sides 



e 



5 



Radius (arcsec) 



10 



15 



Fig. 6. ESO 338-IG04 continued: Upper: Rotation curve (RC) of main component using the "masked" velocity field, 
i.e. with points north-east of the centre and in the eastern tail (up to where the isophotes of the bright central starburst 
region starts in Fig. ||) excluded, and limiting the half sector to S — 30°. As before incl = 55° and PA = 60°. 
Bottom: Rotation curve of the "perpendicular" component, based on S — 20°, incl = 62° and PA = 336° . Its 
centre lies 5" east of the point where the Ha flux has its peak. For further explanation of what is shown in the RCs 
see the caption of Fig. [l] and Sect. 5. 
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Fig. 7. ESO 338-IG04 continued: Based on FP1 data, with lower velocity resolution: Upper left: Continuum isoin- 
tensity image. Note the bright field star superimposed on the western part of the galaxy. Mid-left: Map of the velocity 
dispersion (FWHM) in units of km/s. Upper right: Isovelocity contours of the ionised gas (thick solid lines) superim- 
posed on the monochromatic Ha image (thin solid lines). The spatial scale in each image is shown in the lower right 
corner. North is up, east is to the left. Bottom: Rotation curve (RC) based on the FP1 velocity field and £ = 55°, 
incl — 55° and PA = 60°. The arrows (f) on the abscissa mark the location of 0.4i? 2 5, O.8-R25, I.O-R25 an d l-2i?25- 
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Fig. 8. ESO 338-IG04 B: Upper left: Continuum isointensity image. Mid-left: Monochromatic Ha isointensity map. 
The step between two consecutive levels is the same for the continuum and monochromatic images. Upper right: 
Isovelocity contours of the ionised gas (thick solid lines) superimposed on the monochromatic image (thin solid lines). 
The step between two consecutive contours is 10 times lower than the step used in the mid-left image. The spatial 
scale in each image is shown in the lower right corner. North is up, east is to the left. Bottom: Rotation curve (RC) 
based on S = 40°, incl = 55° and PA = 230°. For further explanation of what is shown in the RC see the caption 
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Fig. 9. ESO 185-IG13: Upper left: Total monochromatic Ha emission (no decomposition of first and secondary- 
components), with the total velocity field overlaid. Mid-left: Profiles within square boxes of 2x2 pixels of the 
decomposed first and second monochromatic components. Upper right: Isovelocity contours of the ionised gas of the 
first component (thick solid lines) superimposed on the first monochromatic Ha image (thin solid lines). The step 
between two consecutive contours is 4 times larger here than those used for the continuum contours. The spatial scale 
in each image is shown in the lower right corner. North is up, east is left. Bottom: Rotation curve (RC) of the primary 
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Fig. 10. ESO 185-IG13 continued: Upper left: Total continuum isointensity image. Mid-left: Line profiles of the 
non-decomposed data plotted within squares of 2x2 pixels. The intensity scale (Y-axis) in each box with line 
profiles is normalised to the peak intensity in that box. Upper right: Isovelocity contours of the ionised gas of the 
second component (thick solid lines) superimposed on the second monochromatic Ha image (thin solid lines). The step 
between two consecutive contours is 4 times larger here than those used for the continuum contours but is the same 
as the step chosen for the first component. The spatial scale in each image is shown in the lower right corner. North is 
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Fig. 11. ESO 400-G43: Upper left: Continuum isointensity image. Mid-left: Monochromatic Ha isointensity map. 
The step between two consecutive levels is the same for the continuum and monochromatic images. Upper right: 
Isovelocity contours of the ionised gas (thick solid lines) superimposed on the monochromatic image (thin solid lines). 
The step between two consecutive contours is 6.7 times larger than the step used in middeft image. Bottom: Rotation 
curve (RC) based on S = 55°, incl = 55° and PA = 225°. For further explanation of what is shown in the RC see 
the caption of Fig. [I] and Sect. 5. 
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Fig. 12. ESO 400-G43 B: Upper left: Continuum isointensity image. Mid-left: Monochromatic Ha isointcnsity 
map. The step between two consecutive levels is 5 times larger for the monochromatic image than for the continuum 
map. Upper right: Isovelocity contours of the ionised gas (thick solid lines) superimposed on the monochromatic 
image (thin solid lines). The step between two consecutive contours is 4 times larger than the step used in mid- left 
image. The spatial scale in each image is shown in the lower right corner. North is up, east is to the left. Bottom: 
Rotation curve (RC) based on5 = 45°, incl = 60° and PA = 305°. The fitted curve (solid line) has been corrected for 
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Fig. 13. Tololo 0341-407: Upper left: Continuum isointensity image. Mid-left: Total Ha emission (addition of both 
components). Upper right: Isovelocity contours of the ionised gas (thick solid lines) of the first component of the 
eastern galaxy, and the isovelocity contours of the second component (thin lines) of the western galaxy superimposed 
on the monochromatic Ha image of the first component (where the eastern galaxy dominates). The step between two 
consecutive contours is the same as those used for the continuum contours. The spatial scale in each image is shown 
in the lower right corner. North is up, east is to the left. Bottom: Rotation curve (RC) of the eastern galaxy first 
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Fig. 14. Tololo 0341-407 continued: Upper left: Individual profiles. Two Gaussian profiles are plotted in each 
pixel, corresponding to the two monochromatic components and the two velocity fields. Mid-left: Total Ha emission 
(addition of both components) with total (non-decomposed) velocity field overlaid. Upper right: Isovelocity contours 
of the ionised gas of the first component of the western galaxy (thick solid lines) and isovelocity contours of the second 
component of the eastern galaxy (thin lines); superimposed on the monochromatic Ha image of the second component 
(where the western galaxy dominates). The step between two consecutive contours is the same as that used for the 



